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Abstract 

A nonstandard systems of differential equations describing two-species phase 
segregation is considered. This system naturally arises in the asymptotic analysis 
recently done by Colli, Gilardi, Krejcf, and Sprekels as the diffusion coefficient in the 
equation governing the evolution of the order parameter tends to zero. In particular, 
a well-posedness result is proved for the limit system. This paper deals with the 
above limit problem in a less general but still very significant framework and provides 
a very simple proof of further regularity for the solution. As a byproduct, a simple 
uniqueness proof is given as well. 
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1 Introduction 

In this paper, we consider the system 

(1 + 2g(p)) d t p + fig'(p) d t p -Ap = (1.1) 

d t p + f\p) = pg\p) (1.2) 

d u p\ T = (1.3) 

/i(0) = po and p(0) = p (1.4) 



*Acknowledgments. Some financial support is coming from the MIUR-PRIN MIUR-PRIN Grant 
2010A2TFX2 "Calculus of variations". 
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in the unknown fields p, and p, where the equations (jl.lj) - (jl.2j) are meant to hold in a 
bounded domain OcK 3 with a smooth boundary T and in some time interval (0, T), and 
where d u in (11 .3p denotes the outward normal derivative. In the recent papers [31 0] , the 
well-posedness of the problem ( ll.ll) - (ll.4l) was investigated, and in particular the existence 
of the solution was proved by considering the system of partial differential equations 
obtained by replacing the ordinary differential equation (jl.2p by the partial differential 
equation 

d t p — aAp + f'(p) = p g\p) with the boundary condition d u p\r = (1.5) 

and performing the asymptotic analysis as a tends to zero. This modified system orig- 
inates from the model introduced in [H], which describes the phase segregation of two 
species (atoms and vacancies, say) on a lattice in presence of diffusion and looks like a 
modification of the well-known Cahn-Hilliard equations (see, e.g., [121 [13]). The state vari- 
ables are the order parameter p (volume density of one of the two species) , which of course 
must attain values in the domain of /', and the chemical potential p, which is required to 
be nonnegative for physical reasons. This system has been studied in a series of papers 
and a number of results has been obtained in several directions [21 [51 El [11] . More- 
over, some of these results hold for a more general system involving a nonlinear elliptic 
operator in divergence form in equation (11.11) . in place of the Laplacian (see [21 HI El fTD]). 
In all of the quoted papers, the function / represents a double-well potential. A ther- 
modynamically relevant example is the so-called logarithmic potential defined (up to an 
additive constant) by the formula 

f(p) = d {p logp+ (1 - p) log(l - p))-c 2 p{l - p) for p G (0, 1), (1.6) 

where c\ and C2 are positive constant with C2 > 2ci in order that / actually presents a 
double well. However, the class of the admissible potentials could be much wider and 
includes both the standard double-well potential defined by 

f(p) = \ (P 2 ~ I) 2 for pGM (1.7) 

and potentials whose convex part is just proper and lower semi continuous, and thus 
possibly non-differentiable, in its effective domain. In the latter case, the monotone part 
of /' is replaced by a multivalued subdifferential and ( 11. 5B has to be read as a differential 
inclusion. In [3], such a wide class of potentials is considered, so that the existence result 
for system (ll.ll) - (ll.4p obtained there is very general. However, the solution constructed 
in this way may be irregular, in principle. Nevertheless, it is unique and a little more 
regular than expected, at the price that the corresponding proofs are rather complicated. 

The present paper deals just with potentials that see example (11.61) as a prototype, 
but it provides simple proofs of further regularity. As an application, we give an easy 
uniqueness proof. 

Our paper is organized as follows. In the next section, we list our assumptions and 
state problem (ll.ip - (ll.4l) in a precise form. In the last section, we present and prove our 
results. 
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2 Assumptions and notations 

We first introduce precise assumptions on the data for the mathematical problem under 
investigation. We assume Q to be a bounded connected open set in M 3 with smooth 
boundary T (treating lower-dimensional cases would require only minor changes) and let 
T G (0, +oo) stand for a final time. We set 

V:=H\n), H:=L 2 {n), and W := {v G H 2 (Vt) : d v v\x = 0}, (2.1) 

and endow the spaces (12.11) with their standard norms, for which we use a self-explanato- 
ry notation like || • \\y. For simplicity, we use the same notation also for powers of these 
spaces. The symbol ( • , ■ ) denotes the duality product between V*, the dual space of V, 
and V itself. Moreover, for p G [1, +oo], we write || • || p for the usual norm both in L P (Q) 
and in L P (Q), where Q := Q x (0, T). For the nonlinearities and the initial data we assume 
that there exist 

p*,p*GM with p*<p* (2.2) 
in order that the combined conditions listed below hold. 

/, g : [p*,p*] R are C 2 functions (2.3) 

g(r) > and g"(r) < for every r G [p*,p*] (2.4) 

f'(p*)<0<g'(j>J and g'(p*) < < f(p*). (2.5) 

Ha G V n L°°(Q) and p > a.e. in Q (2.6) 

Po G V and p* < po < P* a.e. in Q (2.7) 

Notice that the functions /, g, together with their first derivatives, are bounded and 
Lipschitz continuous. Moreover, we remark that the different assumptions of [5] can be 
fulfilled by splitting / as f\ + j'2 with fi nonnegative and convex, and suitably extending 
fx, /2, and g to an open interval including [p*, p*]. In particular, the logarithmic potential 
(11.61) fits the above requirements with p*, p* G (0, 1) and reasonable choices of g. 

Now, we recall the part that follows from the asymptotic analysis performed in [3] and 
is of interest for the present paper. 

Theorem 2.1. Assume that the assumptions ( I2.2p - (l2.7p hold. Then there exists at least 
one pair (p, p) satisfying 

p G L°°(0, T; H) n L 2 (0, T; V) fl L°°(Q) and p>0 a.e. in Q (2.8) 
pG H 1 (0,T;H)nL oo (0,T;V) and p* < p < p* a.e. in Q (2.9) 
u := (1 + 2(7(p))p G ^(O, T; (2.10) 

and solving the problem 

(d t u{t),v)+ [ Vfi{t)-Vv= [ p{t)g'{p{t))d t p{t)v 
Jn Jn 
for allv eV and a.a. t G (0,T) (2.11) 

d t p + f(p)=pg\ P ) a.e. in Q (2.12) 

w(0) = (l + 2g(p ))po and p(0) = po a.e. m fi. (2-13) 
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We observe that the first regularity level for the time derivative of u obtained in [3] 
is d t u G L 4//3 (0,T; V*), that is, a little better than (I2.10p . However, one easily sees that 
d t u G L 2 (0, T; V*) by comparison in f[2TTTj) . on account of (j^ - p^ and (Q. 

Remark 2.2. We can also consider the stronger form of (12.111) . 

/ (1 + 2g(p(t))) dtn{t) v + [ fi(t) g'(p(t)) d t p(t) v + [ V/x(t) -Vv = (2.14) 
Jn Jn Jn 

for all v G V and for a. a. t G (0, T), and observe that it is equivalent to (12. lip provided 
that one can apply the Leibniz rule to the time derivative d t u. This is the case if d t p 
exists and belongs to L 2 (Q). However, Theorem 12.11 does not ensure such a regularity. 
Moreover, (12. lip also includes the homogenous Neumann boundary condition (II. 3p in a 
generalized sense. 

The aim of this paper is to prove that any solution to problem (I2.11I) - (I2.13I) satisfying 
the very mild regularity (I2.8I) - (I2.10I) is in fact much smoother and, in particular, unique. 

Now, we list a number of tools and notations used throughout the paper. First of all, 
we often use the elementary Young inequality, 

ab < ea 2 + — b 2 for every a, b > and e > 0, (2.15) 

and repeatedly account for the Holder and Sobolev inequalities. The precise form of the 
latter we use is the following: 

VcL q (Q) and ||v|| ff < C|H| y for every v G V and q G [1, 6], (2.16) 

where C depends only on Q. Moreover, the above embedding is compact if q < 6, and 
the compactness inequality 

IMI<z < e||Vf || 2 + C 3)e ||w|| 2 for every v G V, q G [1,6), and e > 0, (2.17) 

holds with a constant C q ^ e depending on Q, q, and e, only. Furthermore, we exploit the 
embeddings 

L°°(0, T; H) n L 2 (0, T; V) C L°°(0, T; H) fl L 2 (0, T; L 6 (fi)) C L 10/3 (Q), (2.18) 
as well as the corresponding inequality 

IMIlio/3 ( q) < C(|M|l-(o,7W) + IMIl 2 (o,t ; vo), (2.19) 

which follow from combining the Sobolev embedding V C L 6 (f2) and the well-known 
interpolation inequality for L p spaces. Again, C depends only on Q. Finally, we set 

Q t := Q x (0, t) for t G (0, T] (2.20) 

and use the same symbol small-case c for different constants, that may only depend on Q, 
the final time T, the nonlinearities / and g, and the properties of the data involved in 
the statements at hand. A notation like c e signals a constant that also depends on the 
parameter e. The reader should keep in mind that the meaning of c and c e may change 
from line to line and even within the same chain of inequalities, whereas those constants 
we need to refer to are always denoted by capital letters, just like C in (I2.16P and in (I2.19p . 
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3 Regularity 

In this section, we prove regularity results for the solution to problem (12. lip — (I2.13p under 
the assumption that the conditions (I2.2p — (I2.7D hold (we often avoid writing this). In or- 
der to help the reader, we sketch our strategy. We fix any solution (p, p) to problem 
(I2.1ip -( l2\13p satisfying the regularity requirements ( I2.8p - fl2.10p and recall that all of the 
nonlinear terms involving p are bounded. Moreover, p. is bounded, too (cf. (12. 8p ). Thus, 
(I2.12p implies that even dtp is bounded. Now, we set 

a := 1 + 2g(p) and b \= pg'(p) d t p (3.1) 

and notice that d t a = 2g'(p)d t p. Hence, we have 

a G L°°(0, T; V) D L°°(Q), d t aeL°°(Q), b G L°°(Q), a>l a.e. in Q. (3.2) 

Next, we introduce the associated linear problem 

(d t (az)(t),v) + / Vz{t) Vv= b(t) v for a\\ v e V and a.a. t G (0, T) (3.3) 
Jq Jq 

(az)(0) = (l + 2g(p ))p , (3.4) 

whose unknown z is required to satisfy 

z G L°°(0, T; H) D L 2 (0, T; V) and az G W X,1 {Q,T\ V*), (3.5) 

and observe that z = p is a solution. Then, we prove that (I3.3p - fl3.4l) has a unique 
solution z satisfying ( 13. 5ft . This implies the following. If we regularize ( 13.3p -f l3T4l) and 
perform some a priori estimates on the solution to the regularized problem then these 
estimates still hold for any weak limit. On the other hand, such a limit must be p due 
to uniqueness. This entails further regularity for p. Once the regularity obtained for 
p is sufficiently high, we can even prove uniqueness in a simple way. We observe that 
uniqueness for f l3.3p - fl3.4p is not straightforward. Indeed, (I3.3P is a very weak form (due to 
the very low regularity f)3.5p ) of the homogeneous Neumann boundary value problem for 
the equation 

d t (az) — Az = b , 

which is formally uniformly parabolic. However, the equation is not presented in diver- 
gence form, and a might be discontinuous since no continuity for p is known. At this 
point, we can start with our program. 

Lemma 3.1. Let (p, p) be a solution to f l2.lip -f l2.13p satisfying (I2.8p - (l2.10p . and let a and 
b be defined by (13. ip . Then problem f l3.3p -f l3T4l) has a unique solution z satisfying (13. 5p . 
and this solution coincides with p. 

Proof. Clearly, p satisfies (I3.5P and solves (I3.3p - fl3.4p . As far as uniqueness is concerned, 
we can deal with the correponding homogeneous problem, by linearity. Thus, we fix any 
z satisfying (13. 5p that solves (I3.3p -( l3~4"p . where b is replaced by zero on the right-hand 
side of (13 .3p and the initial value in (13 .4p is zero. Then, we introduce the adjoint problem 
of finding v such that 

v G H\0, T; H) n L 2 (0, T; W) C C°([0, T]; V) (3.6) 
—ad t v — Av = z a.e. in Q and v(T) = 0. (3-7) 
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It is easily seen that (j3.6p - (l3.7j) has a (unique) solution since z G L 2 (0,T; H). We use 
such a solution v as a test function for z, observing that the integration by parts formula 

(d t w(t),<j ) (t))dt=(w(T),<j ) (T))-(w(0),<f ) (0))- f (d t <P(t),w(t))dt 



is actually valid if w G L 2 (0, T; V) fl W lsl (0, T; V*) and G C°([0, T]; V) n tf^O.T; V*) 
Thus, the choice = v is allowed by (13.61) . Hence, we have 

0= f (d t (az)(t),v(t))dt+ [ Vz-Vv 
Jo Jq 

= {{az){T),v(T))-{{az){0),v{0))- f (d t v(t),(az)(t)) dt - [ z Av 



—ad t v — Av)z = / 
' JQ 



JQ 



z 2 



which implies that z — 0. □ 
Corollary 3.2. // z G H\0, T; H) n L°°(0, T; 1/) n L 2 (0, T; so/wes 

a c^z + 9 t a z — Az = b a.e.inQ and z(0) = /x , (3.8) 

i/ien z£ /ioWs z — 

Proof. Indeed, our assumptions on z and (13. 8 j) imply both (13. 5 p and (I3.3p - (l3.4p . Thus, 

the previous lemma gives the result. □ 

Theorem 3.3. Assume that the assumptions (I2.2p -( )2~T1) are fulfilled, and let (fi, p) be a 
solution to (I2.11l) - (l2.13p satisfying (I2.8p - (l2.10p . Then jj also satisfies 

li G H 1 ^, T; H) fl L 2 (0, T; VF). (3.9) 

In particular, n G C°([0, T]; V). 

Proof. Thanks to (13. 2p . there exist two sequences {a n } and {b n } of C 1 functions such that 

a n — >■ a, <9ta n — >■ c^a, fc n — >• b strongly in L 2 {Q) (3.10) 

\a>n\ + \dt a n\ + l&nl < c an d a n > 1 a-e. in Q (3-H) 



Then, we consider for any nGN the following regularization of problem (13. 8p : 



a n d t z n + d t a n z n - Az n = b n a.e.inQ and z n (0) = /i , (3.12) 

complemented with the homogeneous Neumann boundary condition. Due to the regularity 
of the coefficients and assumption (12. 6p on fio, it is easy to see that this problem has a 
(unique) solution z n satisfying 

z n G H\0, T] H) n L°°(0, T; V) n L 2 (0, T; W). (3.13) 
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We add z n to both sides of the equation for convenience. Then, we multiply the resulting 
equality by d t z n and integrate over Q t for any t G (0,T). By owing to (13. lip , we easily 
obtain that 



/ \d t z n \ 2 + ]-\\z n (t)\\y <)- I |V/i | 2 + / (b n + z n - d t a n z n 
JQt 1 1 Jn JQt 

<\ \^^\ 2 + \ \d t z n \ 2 + c[ (l + \z n \ 2 ) 



)dtz r 



<c+t; I \d t z n \ 2 + c I \\z n {s)\\yds. 

1 JQt JO 

By rearranging and applying the Gronwall lemma, we immediately conclude that 

||5fcZn|U 2 (0,T;.ff) + |pre|U°°(0,T;V) < c - 

Moreover, by comparison in (I3.12p . we also find an estimate for ||Az n ||#. So, the above 
bounds and standard regularity results for elliptic equations yield the estimate 

||-Z n ||ffi(0,T;#) + |kn|U°°(0,T;V) + \\z n \\ L 2 (0,T;W) < C. (3-14) 

Therefore, by weak compactness, there exists some z such that 

z n -> z weakly star in H l (0,T; H) n L°°(0,T; V) D L 2 (0,T; W), 

at least for a subsequence. This also implies z n — > z weakly in C°([0,T]; H) and, recall- 
ing (13.1 Op . we infer that 

a d t z + d t a z — Az = b and z(0) = [Aq . 

Now, we apply Corollary 13.21 and conclude that z = //, whence (I3.9P follows. The last 
assertion is a consequence of the well-known embedding ff 1 (0,T;ff) fl L 2 (0,T;W) C 
C°([0,T]-V). □ 

The regularity just established can be improved provided that a stronger assumption 
on the initial datum /i is satisfied, namely 

fi £ W. (3.15) 

We observe that the regularity given in (13.151) implies ( 12. 6 p because of the continuous 
embedding W C C°(Q). The new regularity result is stated in the following theorem, 
whose proof is performed with the same technique as before. 

Theorem 3.4. In addition to (I2.2p - (l2.7p . assume that (13.151) holds, and let (yu, p) be a 
solution to (I2.11l) - (l2.13p satisfying (I2.8p - (I2.10I) . Then [i also satisfies 

dtH G L°°(0, T] H) fl L 2 (0, T; V) and /j, G L°°(0, T; W). (3.16) 

Proof. Theorem 13.31 ensures that (13.91) holds for /x, so that the regularity of d t p, a, and b 
can be updated. Indeed, ( I3.9P implies that (12.121) can be differentiated with respect to 
time and that d 2 p G L 2 (Q). Hence, we also have 

d 2 aeL 2 (Q) and d t b G L 2 (Q), (3.17) 
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which allows us to construct a new approximation of problem ( 13. 8p . Namely, we can 
choose sequences {a n } and {b n } of C 2 functions satisfying A3. lOj) — A3. 11 j) and 

\\<^a n \\ L 2 {Q) + ||<9^„||l 2 (Q) < c for every n. (3.18) 

This leads to a sequence of solutions z n to the corresponding problems A3. 12)) . which, how- 
ever, still keeps all of the properties of the approximation we had established in the proof 
of the previous theorem. But the regularity of the coefficient and the assumption (13. 15)) 
on no also ensure us that z n is smoother and that equation A3. 12)) can be differentiated 
with respect to time. By doing this, we get 

a n d 2 z n + 2d t a n d t z n + d 2 a n z n - Ad t z n = d t b n , 

and we can multiply this equality by d t z n and integrate over Q t , where t G (0, T) is 
arbitrary. By integrating by parts with respect to time the term involving the second 
derivative dfz n , and owing to the inequality a n > 1, we obtain 



I [ \d t z n (t)\ 2 + / \Vd t z n \ 2 
1 Jn JQt 

< \ I a n (0)\d t z n (0)\ 2 - - I d t a n \d t z n \ 2 - [ d 2 a n z n d t z n + [ d t b n d t z n . (3.19) 
1 Jn 1 JQt JQt JQt 

Two of the terms on the right-hand side can be dealt with, by accounting for ( 13. lip . (I3.14p 
and (I3.18p . in the following way: 

- r / d t a n \d t z n \ 2 + f d t b n d t z n <c + c [ \d t z n \ 2 < c. 
1 JQt JQt JQt 

For the first one, we observe that (13 . 1 2j) . (13.151) . and (13.11 1) imply that a n (0) < c and 

\d t z n (0)\ < |6„(0) + Az n (0) - d t a n (0) z n (0)\ < c + | A/x | + c/i <c+\A/j \. 

As /io G W, the integral under investigation is bounded. It remains to handle the third 
term. By using first the Holder and Young inequalities, and then the Sobolev and com- 
pactness inequalities (12. 16)) and (12. 1 T)) with q = 4, we have 



t 

2 



d t z n < / \\d t a n (s)\\ 2 \\z n (s)\\ 4 \\d t z n (s)\\ 4 ds 

Qt 

/ Wn(s)||^+ \\d 2 a n (s)\\ 2 2 \\z n (s)\\lds 



< 

- 2 



1 [ \\Vd t z n (s)\\ 2 H ds + c [ \\d t z n (s)\\ 2 H ds + c [ \\d 2 a n {s)\\l \\z n (s)\\ 2 v ds 

2 Jo Jo Jo 



1 '* 



< o / \\Vd t z n (s)\\ 2 H ds + c, 
1 Jo 

where the last inequality follows from (I3.14p and (I3.18p . By collecting all this and (13. 19)) 
and rearranging, we conclude that 

||eW||z,°°(0,T;iO + \\dtZn\\LHO,T;V) < c. (3.20) 
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Moreover, by comparison in (I3.12p . we also infer that {A-z„} is bounded in L°°(0,T; H). 
Elliptic regularity and ( I3.14p then yield the boundedness of {z n } in L°°(0, T\ W). At this 
point, we use weak compactness once more. We obtain, on a subsequence, that 

z n -> z weakly star in W 1 ' 00 ^, T; H) n H^O, T; V) D L°°(0, T; W) 

and conclude that z — p, as before. Hence, (I3.9P holds, and the proof is complete. □ 



Corollary 3.5. Under the assumptions of Theorem \3.4\ every solution (p, p) to problem 
(l2~TTj) - (12331) also satisfies 

^Gl 10/3 (Q). (3.21) 
Proof. It suffices to combine the first of (13.161) and (I2.18p . □ 

The regularity just achieved allows us to give a rather simple uniqueness proof. 

Theorem 3.6. Assume that (I2.2p - ()2.7p and (I3.15P hold. Then the solution (p, p) given 
by Theorem \2.1\ is unique. 

Proof. We adapt the argument used in [5] to the present situation. We pick two solutions 
(pi,Pi), i = 1,2, and set for convenience 

P '= Px ~ P2 , p:= Pi- p 2 , li ■= g(Pi) , 7 := 7i - 72 , Vi ■= d'ipi) , V ■= Vi ~ V2 ■ 

By accounting for Remark 12.21 we write (I2.14p for both solutions, take the difference, and 
multiply by p. At the same time, we write (12.121) for both solutions, add pi to both sides 
for convenience, move all the nonlinear terms to the right-hand side, and multiply the 
difference by dtp. Then, we integrate and sum up. By owing to the identity 

{(1 + 2jx)d t pi + p\rpd t p\ - (1 + 2'j 2 )d t p2 + ^2^2^ip 2 } /i 
= 2 + 2 Ti) M 2 } + 2 7^A*2 P- + P2T)d t p\ p + P2V2d t pp, (3.22) 
and using the boundedness and the Lipschitz continuity of the nonlinearities, we obtain 

\ I IM*)I 2 + / iv^| 2 + / \d tP \ 2 + \ j | P (t)| 2 

z Jn JQt JQt z Jn 

<c[ \d t p 2 \\p\\p\ + c [ (|p||/i| + |9 t p||^| + |atp||p|). (3.23) 

JQt JQt 

Just the first integral on the right-hand side needs some treatement. By the Holder 
inequality, the compactness inequality ( I2.17P with q = 5, and the elementary Young 
inequality, we have for any e > 

\d t p2\\p\\p\ < [ \\d t p 2 (s)\\ l0/3 \\p(s)\\ 2 \\p(s)\\ 5 ds 
Jo 

< f t \M*)\\l**+ /W2(s)|| 2 0/3llP(*)ll2^ 

Jo Jo 
<e f\\Vp{s)\\ld S + c £ f\\p{ S )\\lds+ [\\d t p 2 (s)\\l 0/3 \\p( S )\\lds. 
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As the function s H- || ( s )|| 10/3 belongs to L 5 / 3 (0,T) by Corollary 13. 5 \ and thus 
to L 1 (0, T), the last integral can be controlled by the left-hand side of ( I3.23P via Gronwall- 
Bellman's lemma (see, e.g., [U Lemma A. 4, p. 156]). Hence, it is sufficient to choose e 
small enough and apply this lemma in order to conclude that p = and p = 0. □ 

Remark 3.7. It is clear that the bootstrap procedure used in the above proofs can be 
continued to provide even more regularity for the solution (p, p) to problem (12 . 1 1 j) — f ]2 . 1 3 j) 
under suitable assumptions on the initial data. 
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